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Abstract 

We construct a self-consistent model which describes a black hole from forma- 
tion to evaporation including the back reaction from the Hawking radiation. In 
the case where a null shell collapses, at the beginning the evaporation occurs, but 
it stops eventually, and a horizon and singularity appear. On the other hand, in 
the generic collapse process of a continuously distributed null matter, the black 
hole evaporates completely without forming a macroscopically large horizon nor 
singularity. We also find a stationary solution in the heat bath, which can be 
regarded as a normal thermodynamic object. 
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1 Introduction 



In the analysis of the black hole evaporation, one usually assumes that a horizon is 
formed in a collapse process, and examines the evaporation and entropy in the static 
black hole PQ-[H]. 

In this paper we try to build a self-consistent model which describes both formation 
and evaporation of a black hole including the back reaction from the Hawking radiation^. 
That is, we solve the semi-classical Einstein equation in a self-consistent manner: 

G, u = 8nG(T^), (1.1) 

where (T^) contains the contribution from both the collapsing matter and the Hawking 
radiation. From the solution we can investigate whether a horizon and singularity are 
formed or not. 

We first consider a null shell as the collapsing matter and construct the geometry by 
connecting the inside flat metric and the outside outgoing Vaidya metric on the shell. 
Note that particle creation generally occurs in a time-dependent gravitational potential, 
and especially, the Hawking radiation can appear without a horizon [22]. We invent a 
formula that evaluates the energy flux of such a process. Then we obtain self-consistent 
equations which determine time evolution of the shell and the radiation. The solution 
shows that the radiation stops, the horizon and singularity appear, and the black hole 
remains forever. 

Next we analyze the case where a continuous null matter collapses and discuss the 
mechanism of the Hawking radiation. It has an onion-like internal structure and evapo- 
rates gradually from the outermost part. Then we write down a self-consistent stationary 
solution in the heat bath. It has neither a macroscopically large horizon nor singularity. 

2 Construction of a model 

We first explain the general idea for construction of a geometry which describes a black 
hole from formation to evaporation. Next we propose a simple model. 

Suppose that a gravitational collapse forms a Schwarzschild black hole as in the left 
of FigJTJ If we take time reversal, the existing black hole goes back to the flat spacetime 
as in the center of FigJTJ Then, if we cover the inside of the horizon and the singularity 
by pasting a collapsing matter, we obtain a geometry which describes both the formation 
and evaporation as in the right of FigJTJ 0. Note that whether this picture is realized 
or not depends on the dynamics. Therefore we need to make some model and solve it 
concretely. 

We will consider the following model. When we take a null shell as the collapsing 
matter, the inside spacetime is flat: 

ds 2 = -dU 2 - 2dUdr + r 2 dQ 2 . (2.1) 

1 Note that we mean by "black hole" not one that has an event horizon defined globally as in the 
rigorous sense, but one that is formed in a semi-classical collapse process. Some authors pursued similar 
ideas [Jg]-[2TJ. 

2 In [3T] a similar diagram is discussed. 
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Figure 1: Penrose diagrams corresponding to the general idea for construction of a 
geometry which describes both formation and evaporation of a black hole. The left one 
represents the formation of a Schwarzschild black hole which has the event horizon and 
singularity. The center one corresponds to the evaporation from the Schwarzschild black 
hole to the flat spacetime. The right one describes both the formation and evaporation. 



As a simple model of the outside metric, we take the outgoing Vaidya metric 

ds 2 = - ( 1 - ^] du 2 - Idudr + rW, (2.2) 



where m{u) = ^jgr is the Bondi mass and the only non-zero component of the Einstein 



tensor is 

n.... = _*! 

~2 



a{u ) 



where the null energy condition implies a < 0@. This is the general spherically symmetric 
metric which satisfies = and G av = except for G uu @. 

Note that the coordinate r must be the same in the both side, because it is defined 
as the radius of 2-sphere and there is no room to rewrite r 2 dQ 2 . On the other hand, u 
is related to U as 

a(u) 



dU = -2dr s =1 V4 du, (2.4) 

r s (u) ' 



3 In [21] and [3S], the ingoing Vaidya metric was used to study the evaporation. 

4 These conditions come from the following discussion. At r > a, where a is the Schwarzschild 
radius of the null shell, we can take Ggg = = because most partial waves with I ^> 1 of the 
radiation do not go through their own centrifugal barrier in V/ ~ . Next the incoming flux can 
be neglected there because of the boundary condition that any energy flow does not come from infinity 
except for the shell. Furthermore if we consider only massless fields, we can assume G M M = because 
the Weyl anomaly vanishes approximately in r > o > ( p . At r ~ a, the ingoing flow and Tgg can 
exist with ! > 1, but we assume to neglect them for the simplest model. Therefore we can consider 
the conditions. In this sense the outgoing Vaidya metric represents the outgoing radiation without the 
gray-body factor. 
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Figure 2: One-shell model. The outgoing Vaidya metric and the flat metric are connected 
by the null shell. The null shell r s (u) approaches its own Schwarzschild radius a(u) 
exponentially. 



where r s (u) is the locus of the null shell. This comes from the fact that r s is an ingoing 
null geodesic in the both sides. Thus a simple model is given by connecting the outgoing 
Vaidya metric and the flat metric with the null shell as in Fig [2j We call it one-shell 
model. 

Here we analyze the locus of the null shell r s (u) for a given function a(u). r s (u) is 
determined by the condition (12.41) : 

dr s (u) _ r s (u) - a(u) 
du ~ 2r s {u) ' ^ ' ' 

This equation tells that the shell will approach its own Schwarzschild radius in the time 
scale ~ a if a(u) changes so slowly that the time scale in which a(u) changes significantly, 
||| , is much larger than a, that is, 

|d| < 1 (2.6) 
Then, in the region r s ~ a, we can replace r s in the denominator with a and solve it as 

u 

r s (u) w a(u) - 2a(u)a(u) + Ca(u)e , (2.7) 

where C is a positive constant. Here the term —lad means that as the shell approaches 
to its Schwarzschild radius in the time scale of 2a, the radius reduces by the evaporation. 
(See FigEJ) Therefore the shell cannot catch up with the radius completely as long as 
a < 0, but it approaches to 

r s {u) a{u) — 2a(u)a(u). (2.8) 

Finally we investigate the surface energy- momentum tensor T£ v on the shell. Using 
the Barrabes-Israel null-shell formalism |27], we estimate 

T£ = (-k ■ w)- 1 <5(r) (Mk"k v + Va^) , (2.9) 
M = — V = - ~ hTs , 2 , (2.10) 
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where v = is the four vector of an observer [v 2 = —1), k^ is the ingoing radial null 
vector which is taken as k^d^ = — \^d u —d r in the Vaidya metric (12. 2p and k^d^ = 2d\j — 

r 

d r in the flat space (12. ip . and a^ u is the metric on the 2-sphere [a ^dx^ dx v = r 2 dQ 2 ). 
The fact that V oc — a(u) > implies that the work done by the shell as it contracts is 
transformed to the Hawking radiation. Thus this model is consistent in energetics. 

Time evolution of this model depends on the functions a(u) and r s (u), so we will 
investigate their dynamics in the following sections. 



3 Flux formula 



We will here construct a flux formula J(u) which, at r ^> a, estimates energy flow from 
the black hole: 

dm . , . 

- = -J(u). (3.1) 

In the Heisenberg picture, we use the Eikonal approximation, the point-splitting regu- 
larization and only the s-wave to obtain 



U{uf 2U{u) 



(3.2) 



U(u) 2 3U(u) 

whose form is the same as the Schwarzian derivative. The derivation is given in the 



Appendix A Note that we can also derive the Planck distribution without horizon (see 



Appendix B) 



First we test the formula in the case without back reaction, that is, in the geometry 
obtained by connecting the Schwarzschild metric and the flat space. In this case, from 
( 12 .7p . r s (u) becomes 

r s (u) = a + Cae~^, (3.3) 
where a becomes constant completely. Then the flux is estimated as 

h 1 



J{u) 



71 



rp2 



(3.4) 



967r a 2 

where T# = j^-. This is the same as thermal radiation from a one-dimensional black 
body with the temperature T#. In this sense, the flux formula (I3.2p is consistent with 
the usual result j2] . Note that this result is the same as the Stefan-Boltzmann law except 
for the coefficient. 

From ( 12 .4p . (12. 5p . (13. ip and (13. 2p . we have obtained the self-consistent equations 
which determine the dynamics of the one-layer model, that is, a(u) and r s (u): 

dr s r s {u) — a{u) 
du 



2r s {u) 



^ = -2GJ{u) = 

du 4tt 

where L = y/ Gh is the Planck length. 



r s [u) 
r.Ju) 



2r s (u) 
3r s (u) 



(3.5) 
(3.6) 
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Figure 3: The numerical result of a(u) in the case of the null shell. The initial conditions 
are given by (14. ip with a(0) = 100 and n = 10/9. 



4 Time evolution of a null shell 

Now we consider the collapse of a null shell by using the one-shell model and investigate 
whether it evaporates or not □. The numerical result of (13.51) and (I3.6P is shown in 
Fig. [31 Here we have chosen the initial conditions given by 

TL — 1 

r s (0) = na(0), r s (0) = — , r s (0) = 0, (4.1) 

In 

where n > 1 is a number, and we assume a(0) ^> l p . The shell does not evaporate 
completely, and a horizon and singularity appear. This asymptotic behavior does not 
depend on the detail of the initial conditions. 

We can understand why the radiation stops in the following manner. Let's recall the 
estimation of the Hawking radiation on the geometry with a = (see (I3.3P ). In that 
case, only the term proportional to contributes to the formula. However, now the 
term a appears in r s , (12. 7p . and the exponential factor will damp for large u. Therefore, 
r s becomes a — 2aa asymptotically as in (I2.8p . Because da is at most of order ti/a <C l p , 
we can approximate 

r s (u) « a(u). (4.2) 
Then (13. 5p and (I3.6P can be solved as 

a(u) = a(0) -B ^'e^^' 2 ; (4.4) 
Jd 

where B and D are integration constants, and B is small and positive. From ( 14. 3 p 
u = 00 corresponds to C, = 00, so ( 14. 4 p shows that a(u) will not necessarily vanish as 
u — > 00. 

5 A similar case was studied in a different set up [55] [25] . 
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Figure 4: Collapse of a continuous null matter to the asymptotic region where (I2.8P 
holds for each shell. The vicinity around a point behaves as in the one-shell model in 
its own time coordinate The redshift factor is required to compare physics between 
different shells. 



Thus we have seen that a collapsing null shell with radius r s radiates for a while, but 
it stops and the radius r s almost stays at the Schwarzschild radius a. Then the horizon 
and singularity appear. A single shell does not evaporate completely even if the back 
reaction from the Hawking radiation is taken into account. 

5 Generalization to a continuous null matter and 
the stationary solution 

We discuss the case where a continuous null matter collapses (see Fig. HJ). 
Let's consider a shell. The metric just outside the shell is given by 

ds 2 = - (l - du' 2 - 2du'dr + r 2 dQ 2 , (5.1) 

where a'(vf) is the Schwarzschild radius corresponding to the total energy of the lower 
shells, and v! is the time coordinate for the shell. The locus of the shell r'(u') follows 

dr'ju') _ rV) - a'{u') 
du' 2r'(u') ' 1 ' ' 

and behaves as in the one-shell model: 

r'iv!) « a' - 2a'^- + C'a'e~^ . (5.3) 
du 

Furthermore, for the shell, the flux formula holds, and we have 

£?=-*"' =-§ {u ' u} ' (5 - 4) 

where J' represents the energy flux measured at infinity if the shells outside did not 
exist. Here we have introduced iV degrees of freedom. In the case of the standard 
model, iV ~ 100 in the energy region higher than the weak scale. Note that introducing 
iV corresponds to replacing l p with \fNl p . 
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5.1 Hawking radiation from each shell 

We will show that the Hawking radiation is emitted from each shell, but only shells near 
the outermost one are relevant because of the large redshift. 
First we introduce a notation 

p' = r' — a' , (5.5) 

which represents distance between the shell r' and the Schwarzschild radius a'. Here we 
assume that the last term in (I5.3j) has already damped for each shell, so that we have 

P > = -2a'^ . (5.6) 

Then the energy flux for each shell depends only on its Schwarzschild radius and does 
not have an explicit u'-dependence: 

S7 - -M ■ («) 

From f)5.6p . p' also becomes a function of a': 

p(a') = 2a f (a) . (5.8) 

Now we consider the junction condition of the adjacent shells. By looking at each 
shell from the both side (see Fig. H]), we obtain 

r' -a' r' - a" 

du > = !! ^ du ". (5.9) 



If a — a" = da is small, we get 

dv! r' — a" p' + da' da' 

— = = ^— = 1 + — . 5.10 

du" r'-a' p' p' K J 

By integrating it, we obtain the redshift factor between a' and a" for finite distance: 

du' f f a> da . 

exp / . (5.11) 



du " ' \ J a" P(a) 
From this equation, the quantity 

d , fdU\ 

^*7 log (*?) (5 ' 12) 

can be rewritten as 

d ( r' {u,) da\ _ da' 1 1 
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where ( 15 .6p and (15. Sp have been used. Thus we have found that for each shell 



e = l„ (5.i4) 

which is independent of a concrete form of f{o!). 

On the other hand, by expressing {v! , U} in terms of we can express the energy 
flow in (15.41) as 



J' = — U 2 -2^ . (5.15) 



From (I5.14p and ( I5.15p . we obtain 



Nh ( 1 4 da'\ 
■^ = 7SZ + • (5-16) 



967r \a' 2 a' 2 du 
By substituting (15.41) and (15.161) iteratively, we have 

/(„<) = J™ _ "X + (a-°) „ i^T" (5.17) 
V ; 967ra' 2 11527rV 4 V ; 6ft H ' V ; 

where = A, which would be the temperature measured at infinity if the shells 
outside did not exist. Thus, any shell can emit the Hawking radiation if the shells are 
continuously distributed so that we can use (15. lip . This result does not depend on the 
behaviour of the shells outside the one we are considering. From (15.41) . (15. 6p and ( I5.17p . 
p' is determined as 

Nl 2 

P(a ' ) = 24^ + ° (a ~ 3) - (5 ' 18) 
By considering the outermost shell, we find that the total Hawking radiation is given 

by 

J(u) = ^rT% + 0(a-% (5.19) 

which coincides with the result for the static Schwarzschild geometry (13.41) . Then by 
applying (15. 4p to the outermost shell, we obtain the time evolution of the size of the 
black hole: 

^ = -2GJ(u) = --^ + O(a-) . (5.20) 
We can also show that the energy spectrum of the radiation follows the Planck distri- 



bution (see Appendix B ). Therefore this black hole evaporates completely as is usually 
expected. However, our model describes how it happens more precisely. Actually the 
black hole evaporates gradually from the outermost shell as if one peels off an onion. 

Here we will check that the total radiation (I5.19P is equal to the sum of the radiation 
from each shell. First let's estimate the radiation from the region between a 1 and a" = 
a' — da as depicted in Fig. HI If there were no shells outside this region, the radiation is 
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estimated as 



f(or)-f(o!-da)e-ti-*-™ « f(a') - [f(a') - ^-da] (l - ^ 

df(a') . , ,/ /\ da 
— —da + j(a)- 



da! p{a!) 
da 
2a'' 

where ^fe ^ da is neglected as a higher term, and (15. 8p is used. By using this and the 
redshift factor, the sum of radiation from each layer is estimated as 



f a da' _ da" 

the sum of radiations = / — e a> p<- a '"> 

2a' 

i r 



da'e p( a '> 



o 



g pO) 



2a 

P(a) 
2a 

~ 2a 

= f(a) = the total radiation . (5-21) 

Here the dominant contribution in the integration comes from the outermost thin region 
with a width about p(a) oc a -1 (see (I5.18P ). Although each shell radiates, the outermost 
region gives the main contribution because of the large redshift. 



5.2 The stationary metric 

We consider the case that the black hole is put in the heat bath with the Hawking 
temperature of the outermost shell for long time so that (12. 8p holds for each shell. It 
is not difficult to calculate the metric for this stationary geometry, and we obtain (see 



Appendix C) 



ds 2 = — -^e dt 2 + -frj^dr 2 + r 2 dQ 2 . (5.22) 

Nl 2 

This expression is valid for r < a + and smoothly connected to the Schwarzschild 

Nl 2 

metric at r = a + ttt^-- This metric does not have a horizon. Here t is the time of the 
flat space at infinity, which is related to the time around the origin T as 



Nl 2 -l^a 2 



dT = —?-e ^ dt . (5.23) 



a- 



This means that T is so much redshifted that T is almost frozen from viewpoint of an 
observer at infinity. Note that this geometry has been obtained self-consistently, so the 
classical limit (h — > 0) does not exist. 
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This metric does not have a large curvature compared with l p 2 in the region r ^> 
Nl p if iV is sufficiently large, N > 100: 

R = > (5-24) 

The singularity around the origin r ~ is controllable in the sense that it can be 
removed by introducing a small shell surrounding the origin. For example, suppose a 
small shell with a ~ C\fNl p comes first, and next, it grows to a large size with a ^> l p 
in the heat bath. Then the outside region r > ao is described as the stationary metric 
(I5.22p . while the center shell is the Schwarzschild black hole with the radius ao which 
does not evaporate forever as in the case of the one-shell model. Therefore we have a 
horizon and singularity around the origin, but their size is small. 

Here we make a comment on the Weyl anomaly. The trace of the Einstein tensor is 
given by 

^ = 6(|| + i). (5.25) 

Because classically the energy-momentum tensor of null shells should be traceless, this 
should be identified with the Weyl anomaly. Actually, if we use the formula of the Weyl 
anomaly for N scalar fields [30], we obtain 

8rf < T ""> - 3 i (l| 4) • < 5 - 26 > 

which agrees with (I5.25P up to numerical coefficients. Therefore the self-consistent 
solution obtained by the Eikonal approximation f)5.22p already contains the effect of the 
Weyl anomaly. 



6 Conclusion and Discussion 

We have solved the semi-classical Einstein equation in a self-consistent manner. We have 
built a model which describes a black hole from formation to evaporation including the 
back reaction from the Hawking radiation. We consider null matter collapse and assume 
that the geometry is obtained by connecting the matter region and the outgoing Vaidya 
metric. 

Using the Eikonal approximation, we have found a formula that gives the energy 
flux of the particle creation in a dynamical geometry. Then we have obtained the 
self-consistent equations which determine time evolution of the collapsing matter and 
radiation. 



10 



As the first example, we have analyzed the case where a single shell collapses and 
solved it numerically and analytically. The shell does not evaporate completely, and a 
horizon and singularity appear. This is not a thermodynamic object but a stable one 
in the sense that it cannot be formed nor evaporated adiabatically in a heat bath. 

Next we have discussed the case where a continuous null matter collapses. Then the 
Hawking radiation occurs not only from the surface but also from the inside. However, 
because of the large redshift, the radiation is emitted substantially only from the region 
around the surface. This black hole evaporates as is usually thought. We then have put it 
in a heat bath and found the stationary metric. It dose not have a macroscopically large 
horizon or singularity. By introducing a small shell around the origin, this singularity 
can be controlled. The metric automatically takes into account the effect of the Weyl 
anomaly. 

There remain some open problems. Our stationary solution has neither horizon 
nor singularity, so the information inside the hole must come back after evaporation. 
However, we don't understand the mechanism clearly yet. For example, suppose that 
we throw a newspaper into the stationary black hole described by our metric. It will 
behave like another null shell going to the hole as it approaches the surface. Clearly 
its energy will be transformed into the Hawking radiation by our mechanism. However, 
the radiation itself comes from the quantum field on the past infinity, or the vacuum. 
How will the information of the newspaper come back? A clue to this problem is that 
we have taken the expectation value of the energy-momentum tensor (T^) in our self- 
consistent equations, which might correspond to the coarse-graining procedure in the 
ordinary statistical mechanics. 

On the other hand, if we put our black hole in a heat bath with the temperature equal 
to the Hawking temperature of the outermost shell, it is completely stationary. In this 
sense, our black hole can be regarded as a thermodynamic object having this temperature 
and its entropy is given by the area law. We don't claim that the information problem is 
solved, but our black hole does not have a macroscopically large horizon and singularity. 
The small singularity around the origin would be resolved by string theory. If it is the 
case, the system is completely well-defined. 
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Appendix A Derivation of the flux formula ( 13.21) 

We will here derive the flux formula (13. 2p by taking only the s-wave and using the Eikonal 
approximation. From (II. ip and (I2.3p . we estimate (T uu ) at r ^> a in the one-shell model 
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(see Fig. [2]). 

First we investigate the behavior of a massless scalar field at r ^> a in the Schwarzschild 
metric: 

ds 2 = -f(r)dt 2 + -^--dr 2 + rW, (A.l) 
fi r ) 

where f(r) — 1 — -. The action of the field <p on this metric is 

s = -J d*xV=j (IfTdptpdM 

= J dtdr *n,m) (dt - dl + V t {r)) <p {l>m) , (A.2) 

l,m 

where we have decomposed the field into partial waves 

^,r,n)=x: yft f' r) M"). 

l,m 

and introduced the new coordinate dr* = y and the effective potential for each partial 
wave with angular momentum I as 

, , p . . f Ul + l) d r f(r)\ 1(1 + 1) 
Vi(r) = f(r) + J^l j „ A_J at r » a . (A.3) 

This implies that only the s-wave survives at r ^> a because partial waves with I > 
have to tunnel their own centrifugal barrier with the rate Pi ~ e~ l . 

Then let's consider the wave equation for scalar field ip on the Vaidya metric (12. 2ft : 

= VV0*0 

-2<9 M + <9 r ^ <9 r y? + - -0 U + ^<9 r U - —ip, (A.4) 

'y J 7" \ IT J 7" 

where I 2 is the Laplacian for angular directions. Here we take only the s-wave 

^r^jw— , (A.5) 
r 

and use the Eikonal approximation (h — )■ 0). Then we get 

= (LZ±d r ^ - 2d u i/?J d r ip . (A.6) 

Therefore, for the outgoing modes, we obtain the equation: 

d r i\) = 0. (A.7) 

Next in this approximation we consider time evolution of the Heisenberg operator 
at r ^> a in the collapsing spacetime. (See Fig. |5j) Before the collapse, the field is given 
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Figure 5: The time evolution of the Heisenberg operator <\> before and after the collapse 
in the Eikonal approximation. 



by the spherical waves: 

f°° doo ( e~ ibJU e +iuju \ 

Mu, r)= — 7= — a " + Ta — a " ' ( A - 8 ) 
Jo Zn yy/Anur y/Airur J 

which corresponds to the field on the flat space. Here the vacuum is defined as the 
Minkowski vacuum: 

a w |0) =0fora;>0. (A.9) 
After the collapse, the field becomes, from (IA.7j) . 

,, N f°° doo ( e - iu > f ^ e +^/H \ 
J Zn \y/A7iujr y/Airur J 

where f(u) is any increasing function of u. We are here using the Eikonal approximation, 
so the phase remains constant on the outgoing mode: 

f(u) = U(u) = -2r», (A.ll) 

where U is the time coordinate in the flat space inside the shell, and at the second 
equality we have used the junction condition on the locus of the shell (12.41) . 

Let's estimate the flux based on the above analysis. We use the point-splitting 
regularization technique to subtract the divergence [3D] : 



(0| : T uu (u) : |0> = lim [(0| : d^d^u') : |0> - (0| : d u( t> Q (u)d u cj> {u') : |0>] , (A.12) 

u— >u 

where the time u is after the collapse, and the r-dependence is not explicitly written 
because r ^> a. First we introduce 

u' = u + e, u" = u - e , (A. 13) 

and expand the equation with respect to e. By using (lA.lOp . the first term in flA.12|) is 
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estimated as 



4vrr 2 (0| : d u <t>{u')d u <t>{u") : |0) 



h 



duo 
2^ 



"OO 7 

ft / f + /V + fV 2 /9. -I 117 — 4- fV 2 /9. + • • • 1p-M2A+ / e 3 /3+- 



— w[/ + /e + fe 2 /2 + •••][/- /e + /e 2 /2 



2ir 

h h 
+ 



3vre 2 8vr 



Z!_2/ 

/ 2 3/ 



+ 0(e 2 



where / = /(m). In the same way, the second term in (1A.12|) is estimated as 

h 



47rr 2 (0| : d u 4> G {u')d u 4> Q {u") : |0) 



87TC 2 



Thus, by using ( lA.llj) . we obtain the flux formula for J(u) = 47rr 2 (0| : T uu (u) : |0) as 



U{u) 2 _ 2U{u) 
U(u) 2 3U(u) 



(A. 14) 



which is (J32]). 



Appendix B Derivation of the Planck distribution 

without horizon 

We emphasize that the Planck distribution can be obtained even if the geometry has 
no horizon. All that is necessary is that the affine parameters on the null generators of 
past and future null infinity are related exponentially [22] . 

In this appendix, we will show that in our model, the expectation value of the number 
of the particle creation takes the form of the Planck distribution with the Hawking 
temperature Th(u) = ^a{u) ' m wmcn a ( u ) changes so slowly that (12. 6p holds. 

We start with reviewing the standard calculation of the Hawking radiation. We 
consider the state in the Heisenberg picture that is annihilated by the positive frequency 
operators in the past infinity a w : 



a w |0) = 0, u > 0. 



(B.l) 



Because the profile of the wave is modified by the gravitational potential, the posi- 
tive frequency operators in the future infinity is a superposition of the positive and 
negative frequency operators in the past infinity a w ', 

bu, = J du'Auurdu' = J du' (Aw'Ow' + ■ (B.2) 
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Then in the future infinity the number operator takes the non-trivial value 

POO 

<0|^M0>= / dJ \A^,\ 2 . (B.3) 



The coefficient A u - U ' is given by the Klein-Gordon inner product: 

A« t - u , = (<p b (u;u),<p* a (U',u;% (B.4) 

where ifb{u; u) is the wave function of the outgoing mode on the future null infinity and 
(p a (U; u) is that on the past null infinity. 



As in |Appendix A[ we will use the Eikonal approximation for the s-wave. Then, 



from ( EOP and (jATTO]) . 



■y £-iuu -y ^,—iuiU 

ip b (u;u)~ , <p a (U;uj) -== . (B.5) 



Then flB.4|) becomes 



CO 



= ^—J — f I due-* M e-* J ' v M. (B.6) 



2tx V w 

Here we need the relation between u and U. 

(1) In the case of a single shell, we can use U(u) = —2r s (u) and r s (u) ~ a(it*) + 

Ca(u*)e 2a (^») where u* is a time when the exponential factor remains. Here note that 
C is positive because r s > a. Then we obtain 

u — u* 

U{u) = -2a(u*) - 2Ca{u*)e~^M. (B.7) 
Thus (IB.6P can be evaluated as 

1 / POO _ u — 11* 

A^ Ut = ^J-l due"*"^^ 5(5I I (B.8) 

where the irrelevant phase factor is dropped. Here the contribution from u away from 
is negligible j2j [22] because the only interval [u* — ka*,u* + ka*\ contributes to the 
integral, where A; is a constant ~ 1. After performing the w-integration, we obtain 

Aj,-u';u, = ^ X [^e-™^"r(i2a(uM, (B.9) 

where the irrelevant phase factor is omitted. 

(2) In the case of the asymptotic region (I2.8p of the continuous matter, we first 
expand a(u) around which is a time in the region: 

{11 — 11 ) 

a(u) = a{u*) + a(u*)(u - u*) = o(w#) — (B.10) 

a[u*y 
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where for simplicity we have normalized the Hawking radiation as a = — 4f- Then the 
redshift factor is estimated as 

dU C a ( u ) da' 1 / \2 1 I \2 i «— tt* 

_ = e ~ Jo ~) = e~^ a{u) = e-^W"*) e + M^). (B.ll) 
du 

Thus we obtain 

U{u) = const. + 2a(w*)e" a ("*) e + ^) = D* + 2C'a(u*)e + ^M , (B.12) 
where C > 0, and have 



•-co,— w';u* n \l / I r, 1/ / 

Z7T V W V-oo ^7T V W 



(B.13) 

This is different from (1B.8|) in the sign of the exponentials, but this integral leads to 



almost the same result: 

A^ I]U , = ^ t £-<««)T(-i2a(«,) W ), (B.14) 
For the both cases, using the formula 

x sinh(7rx) ' 



T{ix)\ 2 = — (B.15) 



we obtain 

2 a{u*) 1 . . 

\A UJ - UJ i ]Ut \ = ; — —g— , (B.16) 

TtLO e T(u t ) _ 1 

where 

TK) = — (B.17) 

By considering a wave packet around «*, we arrive at the Planck-distributed Hawking 
radiation with temperature T(tt*): 

(0|6^|0) u . = — ^ • (B.18) 



Appendix C Derivation of the stationary metric (15.221 ) 



We will derive the stationary metric (I5.22p . The metric ( 15.1]) represents a vicinity around 
a point (V, r') just outside a shell, r' is so close to a' that, from (I5.18p . 

Nl 2 Nl 2 

p> = r i- a i = L^^Ll!iL. (c.i) 

From (I5.4p . ( 15. 6 p and (I5.17p . p'" 1 is also estimated as 

1 2W + 4- (C2) 



p(a') Nl 2 r' 
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By using these, the time coordinate u' is related to that around the origin U as 

du' ( [ a ' da \ r' 2 j^r' 2 . . 

^^{Uw))^;"" • (a3) 

From these, (15. ip can be rewritten as 

r \dU J du 
~ ~ 9A M12 e7HTp dJj2 ~ 2 l^n e ^ P dUdr + r 2 dtt 2 , (C.4) 

where we have replaced r' with r. Here we introduce the time coordinate around the 
origin: 

dT = dU + 24vre~^t r dr, (C.5) 

and express the metric as 

r 2 94-7rr 2 

^"24^ ^ + ^ + (C.6) 
Now we connect it to the outside metric, that is, the Schwarzschild metric at the 

Nl 2 

outermost shell r = a H -• 



24™" 



rfs 2 = -h-^W+fl-^) ~dr 2 + r 2 dQ 2 



-l 



r p 

M * d^ + r^^ + rW. (C.7) 



24vra 2 iV7 2 

Comparing this with (1C.6P at r = a + ~ a, we obtain the relation 

Nl 2 -^a 2 

a 

Therefore we reach the metric 



dT = ~—fe ™T dt. (C.8) 



^ = -24^ 6 " ^ + + rW " (C - 9) 
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